evals %>%
rename(sex = gender,
minority_status = ethnicity) %>%
group_by(sex, minority_status) %>%
summarize(mean_eval = mean(score))# A tibble: 4 × 3
# Groups: sex [2]
sex minority_status mean_eval
<fct> <fct> <dbl>
1 female minority 3.93
2 female not minority 4.13
3 male minority 4.25
4 male not minority 4.23